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Abstract. We have analyzed the trajectories of moving stellar groups in the solar neighborhood in an attempt 

oo 
(N 



to estimate the number of supernova explosions in our local environment during the past 20 million years. Using 
Hipparcos stellar distances and the results of kinematical analyses by Asiain et al. (1999a) on the Pleiades moving 
groups, we are able to show that subgroup Bl, consisting of early type B stars up to 10 M0, but lacking more 
1^-^ . massive objects, has passed through the local interstellar medium within less than 100 pc. Comparing the stellar 

' content of Bl with the initial mass function derived from the analysis of galactic OB associations, we estimate the 

^SJ , number of supernova explosions and find that about 20 supernovae must have occurred during the past ~ 10 — 20 

■ million years, which is suggested to be the age of the Local Bubble; the age of the star cluster is about ~ 20 — 30 

million years. For the first time, this provides strong evidence that the Local Bubble must have been created and 
shaped by multi-supernova explosions and presumably been reheated more than 1 million years ago, consistent 
with recent findings of an excess of ''"Fe in a deep ocean ferromanganese crust. Calculating similarity solutions of 
O ' an expanding superbubble for time-dependent energy input, we show that the number of explosions is sufficient 

, to explain the size of the Local Bubble. The present energy input rate is about Esn ~ 5 x 10^^ erg/s, in good 

^ ' agreement with the estimated local soft X-ray photon output rate. It seems plausible that the origin of the Local 

Bubble is also linked to the formation of the Gould Belt, which originated about 30-60 Myrs ago. 
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1. Introduction Although it has often been suggested that the LB is 

r 1 TT 1 result of one or more supernova (SN) explosions (Cox 

Our solar system IS embedded m a medium of low Hi col- c a j inoo t p xj j. • j- ino,f\ -ti, 

, , „ , , , , nii , , Anderson 1982, Innes & Hartquist 1984), there is no 

umn density (e.g. Frisch & York 1983) largely failed with j- . • > r 1 j-i, • tt n ■ ^-u • ui, 

^ . . . r ~ir f flirect evidcncc lor this hypothcsis. Unlike m the neighbor- 
hot plasma radiating in the soft X-rays (McCammon et . y t u uui -ii • ^ r 

T s , T 1 iiig Loop 1 superbubble, with still ongoing star lormation 

al. 1983). The obvious anticorrelation between neutral gas j u j. ci\t 1 • j. • 4.1, r 4. r-c^ 

, ^. ' ... , , , , and about 40 SN explosions to occur m the future (rLgger 

and A-ray emitting plasma on large angular scales has 1nno^ u • j.- 1 i. r 1 i. ^ -j 

, . ■' ,. r 11111 • 1998), no obvious active cluster of early type stars resides 

been interpreted m terms of an elongated local hot cavity ■ -j .t t-d m 1 • r i- -c • uj. 

„ r. 1 1 r^^^ o 1 mside the LB. Moreover, the analysis of the fane of sight 

( fanaka & Bleeker 1977, Sanders et al., 1977; Snowden et • ii, j- i- r j-u -t a riiv/r i,- i, -t j 

V ' ' ' m the direction of the star p CMa, which extends over 

al., 1990) with an extension between 80 and 200 pc, now r,r,n ^ ■ j- -t j-u j- j-t, -j- • 4. i, 

' ; 1 -r. 1 1 1 / \ 1 200 pc, seems to indicate that the cavity is not homoge- 

commonly known as the Local Bubble (LB). Recently, , nii 1 vi, v -j-*-- 1 ^ri 4. ^ inoc 

r 11 T neously failed with X-ray emitting plasma (Gry et al. 1985, 

based on the spatial distribution 01 stellar EUV sources ttt 1 i , , 1nno^ j i- -ui i- r-tT, ttd 

, , , 1 1 Welsh et al. 1998) and a non-negligible part of the LB ma- 
detected with EUVE and subsequent Nai absorption fane , • , 4- ■ 4- • i, -4-1, 4- 4- c 

^ '- terial must exist m a warm phase with a temperature 01 

studies (Sfeir et al. 1999), Welsh et al. (1999) suggested 1 in4 rri,- 1 4- ttd r 4-- u • 1 

^ 11 1 • 1-1 „ ■ 1 , r- several 10 K. i his rules out a LB formation by one single 

that the LB shows a chimney like structure with no den- oat 1 • • t j • 4. 

SIN explosion since an expanding pressure driven remnant 
nite border towards high Galactic latitude in the northern , en 4.1 -4. • -j 4. 1, u 

would always fail the cavity mside a swept-up shell. 

direction rather than a full bubble. Despite the detailed 

mappings of the LB and our knowledge about its morphol- ^^"^ ^^^'i' ^^^"^^ ^1^° ^^PP^^'^ ^""^ 

ogy, its origin is still a mystery. isotropic nature of the LB emission, coming from EUVE 

observations of cloud shadows in the local ISM. Based on 

Send offprint requests to: T. W. Berghofer these observations, Berghofer et al. (1998) found evidence 
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that the pressure of the hot interstellar medium is the 
same in three different directions near the Galactic plane. 
Therefore, it seems appropriate to distinguish between the 
local cavity as the region largely devoid of Hi gas and the 
present day LB as the X-ray emitting region. These am- 
biguous observational results have consequently led to a 
variety of different ideas concerning the origin of the LB. 

Some time ago the mysterious 7-ray source Geminga 
seemed to be the most promising candidate for a re- 
cent nearby SN explosion, after having been identified 
with an X-ray millisecond pulsar (Gehrels & Chen, 1993). 
However, HST observations have provided a parallactic 
distance of w 160 pc (Caraveo et al., 1996). Therefore, 
due to its high proper motion, Geminga most probably 
originated 3.4 x 10^ yr ago in the Orion region and is thus 
not related to the present day EUV and soft X-ray emis- 
sion of the LB. 

In another scenario, the superbubble picture of the LB 
is disputed, and the LB is conceived as being related to 
a local interarm region between the Sagittarius and the 
Perseus spiral arms of the Galaxy (Bochkarev 1987; Frisch 
1995, 1996). It is assumed that during different epochs of 
star formation in the Sco-Cen association hot gas was gen- 
erated and subsequently swept into the surrounding inho- 
mogeneous medium. While some morphological features of 
the local ISM can be explained by such a model, it cannot 
account for the distinct X-ray shadow observed towards 
the Loop I bubble in the ROSAT PSPC (Position Sensitive 
Proportional Counter) data (Egger & Aschenbach 1995), 
which is most naturally explained by an interaction ring 
between Loop I and the LB. The existence of a "wall" be- 
tween these two bubbles has also been demonstrated by 
stellar absorption line measurements of stars with known 
Hipparcos distances. The formation of such a wall is rather 
difficult to explain once hot plasma is flowing out from 
Loop I directly into the LB. For a detailed discussion we 
refer the reader to Breitschwerdt et al. (2000). 

Recently, Smith & Cox (2001) have performed 1-D 
hydrodynamic simulations to explore whether a low en- 
ergy input rate in combination with a high ambient pres- 
sure can result in the formation of a bubble like the LB. 
They find that 2 — 3 SN explosions within a warm ionized 
medium are sufhcient to explain most observed charac- 
teristics of the LB. It is argued that due to the absence 
of a parent star cluster, these should be due to random 
explosions, e.g. due to runaway OB stars, in the Galactic 
disk. As a result of the low energy input rate and high 
ambient pressure all models predict a collapse of interstel- 
lar bubbles once the evolution time scale is larger than 
about 6 Myr. Although an interesting idea, there is to our 
knowledge as yet no observational evidence of a bubble in 
contraction. 

Given these widely different concepts for the origin of 

the LB, it seems appropriate to investigate the evolution 
of the stellar content of the local ISM back in time for 
the last tens of million years. Thus we are looking for a 
"smoking gun" of a stellar cluster consisting of massive 
stars that may have passed through the LB; its most mas- 



sive members should have already ended as SNe, but the 
later-type stars should still be present. Applying an ini- 
tial mass function (IMF) appropriate for OB stars we are 
able to infer the number of members that may have ex- 
ploded within a region that now forms the local cavity. 
However, these stars may have dispersed and are not eas- 
ily identifiable as a cluster any more. As pointed out by 
Eggen (1996), small perturbations in the space motion of 
stellar clusters can lead to a significant spread of clus- 
ter members on a relatively short time scale. In Sect. 2 
we take a detailed look at the population of stars in the 
solar neighborhood. In Sect. 3, we provide results of our 
investigations on the youngest stars known to exist, and 
in Sect. 4 we derive the time-dependent energy input due 
to successive SN explosions and calculate the evolution of 
the LB. Finally, in Sect. 5, we discuss our findings and 
present our conclusions on the formation of the LB. 

2. The local stellar population 

Investigations of the stellar content in the solar neighbor- 
hood have shown that young stars do exist. Obviously, 
these stars cannot have the same origin as the much older 
stellar population including our Sun. Ongoing star for- 
mation in combination with phase-mixing and relaxation 
must have led to young and old stars now co-existing in 
the solar neighborhood. 

Detailed studies of the Hipparcos data have provided 
the existence of a number of stellar moving groups. Such 
moving groups are star streams identified by a similar ve- 
locity component in the direction of Galactic rotation. A 
complete review of the concept of such star streams and 
the Galactic structure can be found in Eggen (1996). It is 
suggested that some, if not all, of these moving groups 
belong to so-called superclusters. According to Eggen 
(1996), internal and external stresses (e.g., disk heating 
and Galactic differential rotation) may cause perturba- 
tions in the space velocities of the individual stars in a 
cluster, which led the original cluster to expand in a rela- 
tively short time and form a supercluster. 

Recently, Asiain et al. (1999a) studied the kinemati- 
cal properties of 1924 B, A and F main sequence stars 
extracted from the Hipparcos catalog. They included ra- 
dial velocity measurements to derive galactic space ve- 
locity components (U,V,W) and Stronigren photometry 
to determine stellar ages for their sample stars. Using a 
method based on non-parametric density estimators they 
find evidence for moving groups in the four-dimensional 
(U,V,W,age) space. In addition to already known moving 
groups, e.g., the Sirius and the Pleiades moving groups, 
they detect finer substructures with kinematical prop- 
erties apparently consistent with those of nearby open 
clusters or associations. In particular, they find that the 
Pleiades moving group splits into 4 significantly different 
subgroups of vastly different ages ranging from 20 to 300 
Myrs. According to Asiain et al. (1999a), the youngest (Bl 
in their notation) is also the most distinct subgroup and 
the only one still concentrated in space. Most of the 27 
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Fig. 1. Mass histogram of the stars in the Pleiades Bl 
moving group (dark shaded area). The thick solid line 
shows the best fit initial mass function. Since the num- 
ber of stars known so far to be associated with subgroup 
Bl is actually a lower limit, we also plot the uncertainties 
as a hatched region (see Sect. 3). 



member B stars lie in the nearest direction of the Gould 
Belt, I w 300°. The center of mass of these stars is located 
135 ± 15 pc away from the Sun, whereas the distance of 
the Gould Belt in this direction is ^ 160 pc. Obviously, 
this Pleiades moving group is located in the foreground 
of the Gould Belt. According to Asiain et al. (1999b) the 
average age of these stars is 20 ± 10 Myrs. Based on simu- 
lations of the evolution of a stellar complex including the 
influence of the Galactic gravitational potential and the 
effect of disk heating with a constant diffusion coefficient, 
Asiain et al. (1999b) are also able to explain the gross 
properties of the kinematics of the Pleiades moving group 
substructures. 



3. The youngest stars in the solar neighborhood 

The stars associated with Pleiades subgroup Bl consti- 
tute the youngest stars known to exist in the solar neigh- 
borhood. In order to explore whether more massive stars 
existed in this association that already exploded as SNe 
while moving through the LB region, we assume coeval 
star formation and investigate the mass spectrum of the 
stars in the Pleiades subgroup Bl. 

Employing evolutionary tracks by Schaller et al. (1992) 
for the 27 stars in subgroup Bl we derived individual stel- 
lar masses from their spectral types and luminosities. In 
Fig. |l| we show the mass spectrum of the stars in this 
subgroup. As can be seen, subgroup Bl is dominated by 
early B stars of spectral types B3 - B1.5, with masses in 



the range 6-9 Mq. It is worthwhile to mention that 11 of 
these stars have already evolved from the main-sequence. 
Note that the sample of stars investigated by Asiain et 
al. (1999a) is not a complete sample. The apparent deficit 
of stars with masses less than 6 Mq in Fig. |l] can be at- 
tributed to this incompleteness. In fact, the number of 
stars known so far to be associated with subgroup Bl is a 
lower limit. In order to investigate how many stars more 
massive than 9 Mq may have existed in the Pleiades sub- 
group Bl we assume that all stars in the mass range 8-9 
Mq are detected. Next, an initial mass function (IMF) is 
fit to the distribution of stars in the mass range 8-9 Mq . 
We adopted a standard power law IMF with a power law 
index F of the form 



r 



dlog^(logrn) 
dlogm 



(1) 



^ denotes the number of stars per unit logarithmic mass 
interval per unit area. This translates into a number N{m) 
of stars in the mass interval (m, m + dm), 
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N{m)dm = iVo ( — dm. 
\moJ 



(2) 



According to Massey et al. (1995), F = -1.1 ± 0.1 for 
stars in Galactic OB associations with masses in excess 
of 7 Mq. Assuming this value and normalizing the IMF 
to solar masses and the mass spectrum of subgroup Bl, 
N{mo = 8Mo) = 7 (see Fig. 0), we derive 



N(m) =No[ 

^ ' ' Mq 



r-i 



551.6ti«i. 



(3) 



This IMF is shown as a thick solid line in Fig. |T| Note 
that the number of stars known so far to be associated 
with subgroup Bl is actually a lower limit. Therefore, 
we also plot the uncertainties in positive direction as a 
hatched region. The uncertainties take into account the 
uncertainties in the IMF slope (F = -1.1 -I- 0.1 = -1.0), 
which slightly flattens the IMF, and the positive statis- 
tical number errors in the mass bins 8 and 9 Mq. If we 
assume that such a continuous mass distribution is rep- 
resentative for the discrete number of stellar masses in a 
cluster, and that there are no stars if Nim^;^-^) < 1, we 
obtain from Eq. (||), TOmax — 20.2. Therefore, it is unlikely 
that stars with masses in excess of w 20 Mq have existed 
in subgroup Bl. According to Fig. |^ mass depletion be- 
comes apparent for stars with masses m > m^in = 10 Mq; 
the estimated total number of SN events, iVsN, which will 
have occurred in this mass range, is then derived from 



SN 



No X 



21.2 



rn 

w 







r-1 



dm 



(4) 



+ 15.4 
-4.8 ' 



for r = — l.libO.l. Taking into account that two stars are 
still present with masses of 10-11 Mq, 19 SNe are esti- 
mated to have already occurred. Taking into account that 
the sample of stars known to be associated with subgroup 
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Bl may be incomplete, the derived number of SN events 
may be considered as a lower limit. 

We next compare location and space motion of the 
Pleiades subgroup Bl with respect to the LB. Fig. ^pro- 
vides a sketch of the solar neighborhood in a projection 
seen from the top of the galactic plane; the direction to- 
wards the galactic center is to the right (1 = 0°). The po- 
sitions of a number of nearby open stellar clusters (open 
circles in Fig. ^ are shown together with the center of 
mass position of the Pleiades subgroup Bl (filled circle in 
Fig. |2|). We point out that the Pleiades subgroup Bl is 
localized in a small volume in space. The evolution of the 
space distribution of the single member stars is plotted in 
detail in Asiain et al. (1999b, cf. Fig. 4). The LB is ide- 
alized as a sphere; radius and center have been chosen to 
be consistent with the EUVE measurements presented by 
Berghofer et al. (1998, cf. Fig. 13). Additionally, we plot 
the border of the local cavity as derived by Sfeir et al. 
(1999) from Nai absorption line studies. 

The solid line, ending at the actual position of Bl, 
provides the space motion of the moving group during the 
past 30 Myrs in the epicyclic approximation; tickmarks 
indicate the center of mass position 30, 20, and 13 Myrs 
ago. The path of Bl has been computed by employing 
the equation of motion provided by Asiain et al. (1999; 
Eq. (1)). The lAU standard for the solar motion has been 
adopted to transform the Bl velocity components in the 
local standard of rest (LSR) system, following the assump- 
tion that the system of the LB is bound to the LSR. The 
time of -13 Myrs is estimated to reflect the time when the 
most massive star(s) in Bl (M « 20 Mq) must have ex- 
ploded. At this time the center of mass of Bl is located 
inside the volume of the LB. 

We note in passing that there exists a controversy in 
the interpretation of Hipparcos data, resulting in different 
values for the solar motion. While there is good agree- 
ment on Uq, Dehnen & Binney (1998) obtain a value 
Vq « 5.25 kms~^, by arguing that because there is a 
linear relationship between V and the velocity disper- 
sion cr^ (also known as asymmetric drift), which is ful- 
filled for all stars but the youngest, that the latter ones 
should be excluded from the sample of stars. On the other 
hand, Torra et al. (2000) have analyzed the kinematics of 
young stars, using the Hipparcos data plus photometric 
distances, and by taking into account a distortion of the 
velocity field by the expansion of the Gould Belt, obtain 
Vq w 13.4 kms"^, which is close to the lAU standard 
value of Vq ~ 12.0 kms~^. A more detailed discussion of 
these problems is clearly beyond the scope of this paper. 
Assuming the smaller value for Vq by Dehnen & Binney 
(1998) stretches the Bl trajectory to the upper direction 
in Fig. |2| by a factor of ~2; however, the trajectory still 
passes through the LB as outlined by the Nai observations. 

Taking into account the spread in the actual positions 
of individual Bl member stars and the deviations in space 
velocities (AU, AV) = (4.7,3.3) kms^^, a significant frac- 
tion of former more massive members must have crossed 
the LB region even closer to the center, while a few might 
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Fig. 2. Sketch of the solar neighborhood seen from above 
the galactic plane. The center of mass position of Pleiades 
subgroup Bl is labeled with "Bl". The sohd fine, ending 
at the actual position of Bl, provides the trajectory of the 
moving group during the past 30 Myrs in the epicyclic 
approximation (see Sect. 3); center of mass positions 13, 
20, and 30 Myrs ago are labeled with -13, -20, and -30. 
Approximately 13 Myrs ago the most massive Bl star(s) 
(M « 20 Mq) must have exploded. The local cavity con- 
tours as derived from Nai absorption line studies by Sfeir 
et al. (1999) are shown as thick solid lines (dashed lines 
denote directions of uncertain local cavity borders). As 
can be seen, existing Bl member stars (or at least some of 
them, given their spatial spread) should have crossed the 
region, which now forms the Local Bubble. 

have exploded outside. From the present positions of the 
individual stars in Bl and their space velocities we esti- 
mate that an individual star position may deviate from 
the center of mass 13 Myrs ago by up to ^^100 pc. Even if 
the path of Bl did not run through the very center of the 
LB it should be emphasized (and will also be discussed in 
the next section) that the shape or boundary of the bubble 
are mainly determined by the density and pressure of the 
ambient medium rather than by the location of individual 
SNe with respect to the center of explosion. The reason is 
that a shock wave is weak and propagates fast in a rarefied 
hot medium, but moves rather slowly when it encounters 
high density material, which it has to compress and set 
into motion. 

4. Energy budget and evolution of the Local 
Bubble 

While the details of the dynamics and evolution of the LB 
have to be deferred to a future paper, which concentrates 
on numerical simulations, we will discuss in the following 
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the time-dependent energy input rate and its consequences 
for the evolution of the LB. 

Obviously, the contribution by stellar winds is only sig- 
nificant at early times and thus is relevant only if massive 
O stars exist in the cluster, which is not the case here. We 
will therefore only consider the input rate by SN events, 
which is by far the dominant at later times (Mac Low & 
McCray 1988). The main sequence life times, Tms, of stars 
within the mass range 7 Mq < m < 30 Mq can empiri- 
cally be approximated by t^s = 3 10^ (to/[10 Mo])"" yr 
(Stothers 1972), with a = 1.6. Since this defines m as a 
function of time r, implicitly assuming that the energy 
input can be described as a continuous process, we obtain 



TO 

m;; 



-1/c 



(5) 



with C = 3.762 x 10^^ s. It should be borne in mind that 
describing discrete events (i.e. SN explosions) by a con- 
tinuous mass distribution function in time may disregard 
considerable fluctuations in the time-dependent quantities 
calculated below. However, the results obtained should be 
accurate enough in a time averaged sense. Let Lsb be the 
energy input rate per unit time due to a number of suc- 
cessive SN explosions but with a constant energy input of 
i?SN = 10^^ erg each, then the cumulative number of SNe 
down to mass to reads 



iVsN(TO,) = 501.45 



Mr. 



B. 



(6) 



with B = 18.58, and 



Lsb = -ttC^sn-E'Sn) — Esn — 37— 
at at 



SN- 



dNsjq dm dr 
dm dr dt 



using Eqs. (Q) and 

Lsb = 344.75 C^/"i"(^/"+i) 



(8) 



with the above values for F and a, we obtain the handy 
formula 



^SB 



(9) 



where Lo = 4.085 x lO^^erg/s, 6 = -{T/a + 1) = -0.3125 
and tj — t/10^ yr. Here and in the following, the variable t 
has been used as the dynamical time for energy input into 
the LB and its evolution. Assuming coeval star formation 
we estimate from the lack of stars with to > 10 — 11 M0 
the age of the star cluster to be Td < 2.5 x 10'' yr; good 
agreement with the evolution tracks by Schaller (1992) is 
obtained for Td — 2.3 x 10'' (s. Fig. ||). The present dy- 
namical age for the bubble evolution is therefore given by 
to = Tci - To = 1.3 x 10'' yr, where tq = Tmsirrimax) ~ 9.9 
Myr defines point zero as the main sequence life-time for 
the most massive star {mmax ~ 20 M0). The dynamical 
time is then given hy t — t — To, and hence dr = dt. 

The assumption of a continuous energy input into the 
ISM at the center of the first explosion is much better 



than thought of at first glance. The reason is that every 
off-center blast wave in a hot rarefied medium with a high 
speed of sound rapidly degenerates into a Mach 1 shock 
with a low factor of compressibility of the downstream 
medium. Therefore, hardly any material is swept up and 
the bubble remains roughly uniform with negligible radia- 
tive cooling due to its high initial temperature. The net 
effect is simply pressurization of the bubble. The time for 
the interior weak shock to reach the bubble boundary is of 
the order of a sound crossing time and thus during most 
of the evolution much smaller than the dynamical time 
scale; hence the eccentricity of subsequent explosions has 
only a minor infiuence on the evolution of the bubble. 

Similarity solutions of superbubbles (e.g. McCray & 
Kafatos 1987) and their classical analogue, stellar wind 
bubbles (e.g. Pikelner 1968, Dyson & deVries 1972, Castor 
et al. 1975, Weaver et al. 1977), can be calculated, pro- 
vided that the dependent variables and boundary condi- 
tions do not explicitly depend on time or length scales 
that may enter the problem. This is still true for the case 
of a time-dependent energy input rate, Lsb = Lot^ and a 
density gradient of the ambient medium, p{r) = por^^. 

Similarity conditions hold during the adiabatic expan- 
sion phase as long as the bubble pressure is much larger 
than the ambient medium pressure, i.e. Pb 3> Po- The 
mass of the bubble is dominated by the mass swept into 
a shell, 



Mshir) 



p(r')dV' 



(10) 



and the energy input is shared between kinetic and ther- 
mal energy, 



(7) Ethir) = l/ij-l) I p{r')d\' 



(11) 



Using 7 ~ 5/3 for the ratio of specific heats, observing 
that Pfc remains uniform, and applying spherical symme- 
try, conservation of momentum 



-{MshRb)=^TTRlPb, 
dt 

and energy 
dEth 



dt 



= LsBit) - inRtRbPb , 



yields the solution 
Rb = Af" 

A = 



(5-/3)3(3-/3) 



{7S-I3-S/3 + 11)(4(5 + 7 - (5/3 - 2/3) 



(12) 



(13) 



(14) 

1/(5-/9) 

X 



Lo 



l-i- = 



27r((5 + 3)po 
(5 + 3 



1/(5-/3) 



(15) 
(16) 



since the swept-up shell is usually thin, the bubble and 
shell radius can be treated as equal during the energy 
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driven phase and are denoted by Rt- The similarity vari- 
able in the case considered here is given by 



r/a 



5-/3 



(17) 



To keep matters simple, the ambient density is assumed 
to be constant {j3 = 0), although it is possible that 
some of the explosions occur inside a molecular cloud (s. 
Breitschwerdt & Schmutzler 1994). With the previously 
used values, the radius of the bubble evolves as 



Rb = 251 



2 X 10"24g/pj^3 

Po 



1/5 



.0.5375 



pc. 



(18) 



and the expansion velocity as 



Rb = 13.22 (^^^^^^f^) t^"-'''' km/s . (19) 

As a consequence of a decreasing energy input rate the 
expansion law exponent in Eq. ( p^ ) is between a Sedov 
{ji = 0.4) and a stellar wind (/i = 0.6) type solution. 
According to Eqs. (|l^) and (|l^) the radius of the LB will 
be 289 pc and 158 pc and its velocity is 11.7 km/s and 6.4 
km/s, if the ambient density is po = 2 x 10~^'*g/cm'^ and 
Po = 4 X 10"^'^ g/cm'^, or no = Icm""^ and no = 20cm~^, 
respectively. In the latter case the value of the ambient 
density would correspond roughly to that of the cold neu- 
tral medium. It is interesting to note that the radius, the 
dynamical time or age of the LB, and the rate of SNe 
ifsN ~ 1/(6.5 X 10^ yr)) are very similar to that of the 
neighboring Loop I superbubble; Egger (1998) estimated 
its radius to be Rli = 158 pc, its dynamical age to be 
tdyn ~ 10^ yr and the SN rate to be fsN 1/(6 x 10^ yr). 
Loop I is an active superbubble and has been powered 
by about 40 SNe in the Sco-Cen OB association, which 
is about twice as many as in the case of the LB; in this 
respect, its radius is also fairly small. 

According to ROSAT PSPC broad band spectral fit- 
ting the Loop I bubble temperature is about 2.5 x 10^ K 
(Egger 1998). This gives support to the idea that we live 
in an environment of ongoing star formation, although the 
LB itself is probably extinct. 

There are several reasons why we may have overesti- 
mated the size of the LB in the similarity solutions above. 
Firstly, the mass inside the bubble is significantly higher 
than the pure ejecta mass, as can be inferred from the 
ROSAT X-ray emission measures; when assuming bubble 
parameters of Rb = 100 pc and rif, = 5 x 10"'^ cm~'^ (e.g. 
Snowden et al. 1990) a mass of at least 600 Mq is de- 
rived whereas non-equilibrium ionization plasma models 
(Breitschwerdt & Schmutzler 1994) result in a more than 
a factor of five higher mass. 

The contribution of ejecta is of the order of 100 Mq, 
and the bulk of the bubble mass is due to hydrodynamic 
mixing of shell material, heat conduction between shell 
and bubble and evaporation of entrained clouds; therefore 
the flow is mass-loaded. The net effect is to reduce the 



amount of specific energy per unit mass, because the mate- 
rial mixed in is essentially cold, and thereby also increases 
the rate of radiative cooling. Secondly, the stellar associ- 
ation has probably been surrounded by a molecular cloud 
with a density in excess of uq — 100 cm^'^ with subsequent 
break-out of the bubble and dispersal of the parent cloud 
(Breitschwerdt et al. 1996). Thirdly, the number of SN ex- 
plosions could be less; here we have assumed that all 20 
SNe have occurred inside the LB. This need not be the 
case as is also suggested by Fig. ^. ROSAT PSPC obser- 
vations have revealed an annular shadow centered toward 
the direction (^n = 335°, bu = 0°), which has been inter- 
preted as an interaction between the LB and the neigh- 
boring Loop I superbubble (Egger & Aschenbach 1995). 
The trajectory of the cluster Bl may have partly crossed 
the Loop I region. Alternatively and more likely, part of 
the thermal energy might have been liberated into the 
Galactic halo, since there is some evidence (see Sect. |5|) 
that the LB is open toward the North Galactic Pole. It 
should also be mentioned that due to small number statis- 
tics the true number of SNe can vary by a factor of 2. 
Finally, although there is no stringent evidence, it would 
be very unusual, if the LB would not be bounded by a 
magnetic field, whose tension and pressure forces would 
decrease the size of the LB. 

Given these uncertainties and the fact that the sim- 
ple analytic model discussed above can be considered as 
an upper limit, the agreement with observations is quite 
good. The bubble radius and shell velocity are rather in- 
sensitive to the energy input rate and the ambient den- 
sity (due to the power of 1/5) and therefore not well con- 
strained, but depend more sensitively on the expansion 
time scale. Thus we can assert with some confidence that 
the age of the LB should be between 1 — 2 x 10^ yr. 

The radiative cooling time scale of the LB can be 
estimated using Kahn's (1976) cooling law (in the tem- 



perature range 10^ < T < 10^ K), tc 



,3/2 



I with 



P/p^/^ being the adiabatic parameter and q 



4 X 1032(,m6/(s4g)^ Using n,, = 5 x 10"^ ^^d Tb 



10^ 



K, which have been derived from X-ray broad band spec- 
tral fitting of a plasma in coUisional ionization equilibrium 
(e.g. Snowden et al. 1990), we obtain tc ~ 1.3 x 10^ yr, 
comparable to the dynamical time scale. If the X-ray emis- 
sion is mainly due to the delayed recombination of highly 
ionized species of a plasma not in in ionization equilib- 
rium, then the cooling time can be as short as 6 x 10^ 
yr (Breitschwerdt et al. 1996). In either case, the energy 
input of previous SNe has largely been used up for expan- 
sion and radiative cooling of the bubble. We thus conclude 
that the emergent X-ray emission must be largely due to 
the energy input of the last supernova. 

We can express the number of stars (and potential 
SNe) in the mass interval (m, m+dm), N{m), by a number 
of stars ending their lives in the time interval (r, r + dr), 
N{t), using Eq. (||) and t = t + tq (tq « 10*" yr), which 
results in 



N{t) 



9.69710-20 X (T/[sec])('-^^/" 
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= 1.03 X 



1.3125 



lO^yr, 
L03 X (1+^7)1-3125 



(20) 



The average energy input rate by recent SNe, E'en, (not 
to be confounded with the cumulative rate of Eq. (|^)) is 
then given by 



EsN — E; 



dN{T) dr 



dr dt 

1.2710^2 X (T/[sec])^^-"-^)/" 
4.32 X 10^^(1 + t7)"-3i^^erg/s. 



(21) 



As a result of the increasing number of exploding stars 
with time (compensating their increasing main sequence 
life time) and the assumption that the individual contri- 
bution is constant, i.e. E'en = lO^i erg, the energy input 
rate (using L = —1.1) is a moderately increasing function 
with time. 

Analysis of ferromanganese crust samples in deep 
ocean layers have shown an increased flux rate of ^"Fe, 
which is indicative of at least one SN which exploded in- 
side the LB at about 5 Myr ago (Knie et al. 1999). Thus, 
using tj = 0.8, we obtain for the associated energy input 
rate, Esn = 5.18 x 10^^^ erg/s. This is a moderate value, 
but still comfortably larger by a factor of 5 than the total 
present day soft X-ray photon energy output rate of the 
LB. 



5. Discussion and Conclusions 

The fact that the Gould Belt stars show receding veloci- 
ties, and that the stellar associations are also closely re- 
lated to molecular clouds and local Hi (see Poppel 1997) 
argues in favor of a common history of gas and stars in 
the local environment. The nature of the Gould Belt is 
yet unclear; current ideas favor either an explosive origin 
(Blaauw 1965, Olano 1982) or a collision between a high 
velocity cloud and the Galactic disk (e.g. Tenorio-Tagle et 
al. 1986; Comeron & Torra 1992). It is noteworthy that 
the LB extends perpendicular to the Gould Belt rather 
than to the Galactic disk (Sfeir et al. 1999), hinting at an 
origin of the LB connected to the Gould Belt stars and 
gas. 

Regardless whether the expansion was caused by mul- 
tiple SN events or by the oblique impact of a high velocity 
cloud, the implications are similar: the shock fronts that 
have propagated into the ISM, presumably at some angle 
with respect to the disk, will overrun ambient clouds. If 
the clouds have low masses, they could be disrupted hy- 
drodynamically and/or evaporated, thereby mass loading 
and decelerating the flow. As a matter of fact, the spec- 
tral properties of the soft X-ray emission combined with 
the Hi distribution in the local cavity show that the hot 
LB plasma must consist of a total mass of the order of 
1000 Mq. This cannot be provided solely by ejecta ma- 
terial, but must be due to hydrodynamic mixing and/or 
evaporation of ambient gas clouds. In the case of more 



massive clouds, the shock propagating into the gas will 
compress it without destroying the cloud, so that a grow- 
ing number of Jeans unstable clumps will eventually col- 
lapse into protostellar cores. A considerable number of 
nearby stars with ages substantially lower than the Gould 
Belt's kinematical age and the existence of active super- 
bubbles (e.g. Loop I with an estimated dynamical age of 
10^ years, see Egger 1998) seem to support this hypothesis. 
Note that the age of the Gould Belt is rather uncertain, 
Comeron (1999) provides an age estimate of 34 ± 3 Myrs 
whereas Torra et al. (2000) derive an age in the range of 
30-60 Myrs. 

In the above picture, the origin of the LB is con- 
nected to the evolution of the Gould Belt itself as the 
disturbed medium that formed the LB. We favor a sce- 
nario in which the LB was created by about 10 — 20 
SNe about 1 — 2 x 10'' yr ago, with a preferred value 
of 1.3 X 10^ yr; the star cluster itself is about 2.3 x 10^ 
yr. It is conceivable that the stars have also been embed- 
ded in a molecular cloud for some time (Breitschwerdt 
et al. 1996). The present day total LB X-ray luminos- 
ity is fairly moderate, Lx ~ 10^^ erg/s at most (Cox & 
Reynolds 1987), and can easily be supplied by the last SN 
reheating the cavity about a million years ago. This could 
also explain why only part of the local cavity is illumi- 
nated now by X-rays. Note that the sound crossing time 
out to a distance of 200 pc (in direction of /3CMa) in a 10^ 
K plasma is roughly 2 Myrs. Moreover, the recombination 
time scale of major coolants, such as Ovii and Cvi can be 
estimated from analytical approximations to the rate co- 
efficients (ShuU & van Steenberg 1982) to be 4.7 and 12.2 
Myrs, respectively, for commonly used LB fit parameters 
of nf, = 5 X 10-3 cm-3 and Tf, = 10^ K (Breitschwerdt & 
Freyberg 2002). Therefore the emergent X-ray spectrum 
should be a fairly complex superposition of freshly ion- 
ized and recombining plasma (from previous explosions) 
and definitely be out ionization equilibrium. 

We have shown that the subgroup Bl of the Pleiades 
moving group has actually passed through the LB. If the 
IMF of this cluster of young stars is similar to the OB as- 
sociations we observe in our local neighborhood, then ob- 
viously some of these stars must have exploded as SNe on 
their way. Since the kinematic age of Bl can be fairly well 
determined (assuming they are not bound gravitationally) 
to be TBI ~ 20 Myrs, it can certainly not be responsible 
for the origin of the Gould Belt itself. Note that all age 
determinations for the Gould Belt point to ages older than 
30 Myrs. However, the SN events in Bl should have left 
their imprint on the LB. One important question is when 
the last one injected mass and energy. According to Eq. 
(H) and Fig. |l| about 4-5 stars of 10 - 11 Mq should have 
been members of Bl. Since two of them are still with us 
and the average life time of such stars is about 25 Myrs, 
while stars with 9 - 10 Mq live 30 Myrs the last SN is 
expected to have exploded about 1 — 2 Myrs ago, or 3 
SNe occurred within the last 5 Myrs. This is roughly con- 
sistent with the measurements of radioactive ^°Fe (Knie 
et al. 1999), where the high flux rate is interpreted by 
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one SN, but within a radius of 30 pc from the Sun and a 
progenitor mass of about 15 Mq. 

The existence of a "Local Chimney" as claimed by Sfcir 
et al. (1999) and Welsh et al. (1999) may support the 
hypothesis that the LB is created and shaped by as many 
as 20 SNe, because it allows a significant fraction of the 
large superbubble thermal energy to be vented into the 
Galactic halo. Consequently, after the last explosion, the 
LB is left in a state of long-term recombination. According 
to non-equilibrium models (Breitschwerdt & Schmutzler 
1994; Breitschwerdt 2002) the observed modest amount 
of soft X-ray emission can still be generated, even if the 
last SN occurred more than a million years ago. 

Bl still provides another 10 SN candidate members. 
Since it is about to leave the volume of the local cavity, 
the LB will remain in a state of recombination, which 
is what we probably observe at present. Reheating will 
become progressively less important and the future of the 
LB should be fairly quiescent. 
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